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Abstract. The production of charmonia in heavy-ion collisions is investigated 
within a kinetic theory framework simultaneously accounting for dissociation and 
regeneration processes in both quark-gluon plasma (QGP) and hadron-gas phases 
of the reaction. In-medium modifications of open-charm states (c-quarks, D- 
mesons) and the survival of J/i/i mesons in the QGP are included as inferred from 
lattice QCD. Pertinent consequences on equilibrium charmonium abundances are 
evaluated and found to be especially relevant to explain the measured centrality 
dependence of the ip' /tj) ratio at SPS. Predictions for recent In-In experiments, 
as well as comparisons to current Au-Au data from RHIC, are provided. 



1. Introduction 

Heavy-flavor bound states constitute a valuable probe of the hot/dense strongly 
interacting matter formed in relativistic collisions of heavy nuclei. It was first 
suggested in jl] to identify J/iJj suppression as a signature of a deconfined medium 
due to color Debye screening, with tightly bound cc states presumably being robust 
in a hadron gas (HG). However, more recently it has been realized that c-quark 
reinteractions in the medium can lead to regeneration of charmonium states through 
c-c coalescence [21E]: especially if charm production is abundant {e.g., iVcc~10-20 in 
central Au-Au collisions at RHIC). 

Further insights into charm(onium) properties at finite temperature T have re- 
cently been provided by lattice QCD (LQCD) calculations which indicate (i) a 
continuous reduction of the open-charm threshold with increasing matter tempera- 
ture 0] and (ii) the survival of low- lying charmonia (77c, J/ip) up to ~ 2 Tc |S1E1IZ|- 

We here present an approach to charmonium production at SPS and RHIC [HI in 
which in-medium charm properties are modeled in accord with LQCD results and 
implemented into a kinetic rate equation, solved for a schematic thermal fireball 
expansion ,9,. It enables a simultaneous treatment of charmonium dissociation and 
regeneration throughout the evolution of the system, thus improving on our earlier 
constructed "two-component" model ^3 E] where suppression in QGP and HG was 
combined with statistical production at hadronization. 
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2. In-medium properties of open and hidden charm in equiUbrium 

LQCD calculations at finite T show that the free energy of a static heavy-quark pair 
(QQ), Fqq—V-TS, reaches a (constant) plateau at large spatial separation which 
decreases with increasing T, continuously across the phase transition. We associate 
this behavior with a decreasing open-charm thresholdf; in the hadronic phase, a 
plausible interpretation is provided by the (partial) restoration of chiral symmetry: 
as T increases, we expect a reduction of constituent light-quark masses, inducing 
smaller D-meson masses. The evaluation of the former, within a NJL model at 
finite temperature and density, typically amounts to Am(Tc) ~ —140 MeV. Requiring 
continuity across Tc, the effective charm-quark mass in the QGP is around m* ~ 1.6- 
1.7 GeV, which might be associated with a thermal correlation energy of c-quarks. 

At a fixed number, Ncc, of cc pairs in the system (as expected for production 
exclusively in initial N-N collisions), in- medium modified open-charm states affect the 
(thermal) equilibrium abundance of charmonia given by the densities n^(r, 7c) = 
rf$7c / ^|^/*(m^i,, T), where /* are Bose distribution functions (using vacuum 5* 
masses), and is a spin-isospin degeneracy. The charm-quark fugacity 7c encodes 
chemical off-equilibrium effects being adjusted to Ncs through 

UV ic op; *=,,^,J/^,... 

where A'op = Vnop{'m* jj;T) denotes the total equilibrium number of open-charm 
states in either QGP (c, c quarks) or HG phase (charmed hadrons). In the QGP, an 
in-medium c-quark mass, larger than the perturbative value of nic ~ 1.2 GeV, leads 
to an increase in \1/ equilibrium abundances since it becomes energetically favorable to 
distribute cc pairs into states, especially if m<p < 2m*. Conversely, in the hadronic 
phase, reduced in-medium £)-meson masses lower the ^ equilibrium level. 

3. Rate equations in heavy-ion collisions 

To make contact with the dynamical situation encountered in heavy-ion collisions, we 
employ a kinetic theory framework. If open-charm states are in thermal equilibrium, 
the evolution equation for the number, N^j,, of charmonium states present in the system 
is given by = [-/V^ — N^] . The charmonium equilibrium numbers are 
determined as described above, except, however, for additional thermal off-equilibrium 
corrections expected in heavy-ion collisions, namely: (i) a thermal relaxation time 
which effectively reduces during the early times of the collision (see ^Jj for 
details), and (ii) a correlation volume V^orr implemented into the argument of the 
Bessel functions in Eq. to account for the locality of c-c production. The thermal 
widths of charmonia, F^r = (r^)~^, are obtained by convoluting their inelastic cross 
sections with thermal distributions of the matter constituents. In the QGP phase, we 
use parton-induced "quasifree" breakup reactions ^UIj i + '^^i + c + c {i = g , q, q) , 
accounting for the different binding energies of the respective charmonium states {tp, ip' 
and x)- In the HG, we compute J/V' breakup by pions and rhos within a flavor- 
SU{4:) effective lagrangian formalism (see jj^ ^'^d references therein), augmented 
by geometric scaling for the ^p' and x- Here, the consequences of in-medium D- 
meson masses are two-fold: first, inelastic reaction rates increase due to a larger 



t Note that for a quantitative assessment, the entropy term in Fqq should be removed. 
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Figure 1. Upper left panel: t/j'/i/' ratio in P6(158AGeV)-Pfe at SPS; solid 
(dashed) line: with (without) in-medium effects using a uniform nuclear 
absorption cross section (Tnuc = 4.4 mb; dash-dotted line: as solid line but with 



an updated value for the ip' nuclear absorption cross section, tr, 



7.9±0.6 mb. 



Lower left panel: J/V'/^coll a-t mid-rapidity vs. Afpart in 200 AGeV Au-Au at 
RHIC. The total yield (full curve) is dominated by regenerated J/i/i's (dash-dotted 
curve) with almost no direct J/i/>'s remaining (dashed curve). The band reflects 
the uncertainty on the magnitude of in-medium effects. Upper (lower) right panel: 
centrality dependence of the tp' /ip {JI-ip/DY) ratio in /n(170 AGeV)-/?i at SPS. 



available phase space and second, the DD threshold can move below the charmonium 
mass, which additionally enables direct decays ^' — > DD The rate equations are 
supplemented by initial conditions, N^, for which we take experimental yields from 
pp extrapolated in centrality by N-N collision scaling and "pre-equilibrium" nuclear 
absorption. The rate equations are then integrated over the space-time history of the 
collision according to a schematic thermal fireball evolution 

4. Comparison with experiments 

Pb-Pb at SPS — At SPS energies {y/JNN=17.3 GeV), primordial charmonium 
production is large compared to charmonium equilibrium abundances, 3> N!^"^ , 
implying little regeneration. As shown in |H], the centrality dependence of the J/ip 
over Drell-Yan ratio in Pb-Pb collisions is well reproduced and QGP formation is 
characterized by J/tp suppression. The consequences of in-medium effects at SPS are 
particularly pronounced in the i^'/ip ratio, cf. the upper left panel of Fig. ^ With 
vacuum _D-meson masses (dashed line) our calculation underestimates 4'' suppression. 
The calculation including medium effects (full line) improves the agreement with NA50 
data mi substantially, which is a direct consequence of the reduction of the DD 
threshold in the HG, opening the -0' — *■ DD decay channel. The -0' data set (including 
p-A collisions) has recently been reanalyzed by NA50 QS| (diamonds), deducing a 
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stronger nuclear absorption of the ij]' , (T„uc(V'') = 7.9 ± 0.6 mb > o'nuc('0) = 4.4 nib. 
Our calculation with the correspondingly updated values of tJnuc is shown by the dash- 
dotted line, confirming the need for in-medium effects to reproduce the ip' ratio. 

Predictions for In-In at SPS — The NA60 collaboration will continue to study 
charm(onium) production at SPS and we present in the right panel of Fig. ^ our 
predictions for the V'VV' and J/ip/DY ratios as a function of centrality in /n(170 
AGeV)-/n collisions. Similar to the Pb-Pb system, J/tfj regeneration (dash-dotted 
curve in the lower right plot) is rather limited and the main effect remains J/V' 
suppression (dashed curve) exhibiting a marked departure from nuclear absorption 
(dotted line). Medium effects are still appreciable, as illustrated in the upper right 
panel where a gradual decrease of the tp'/^J ratio is predicted as a function of A^part- 

RHIC results — Our calculations at full RHIC energy {^Jsnn = 200 GeV) are 
compared to published PHENIX data ^H] in the lower left of Fig.^ Contrary to SPS, 
the J/il) yield in central Au-Au collisions (full curve) is dominated by regenerated 
J/tjj's (dash-dotted curve) while primordial J/tjj's are almost completely suppressed 
(dashed line) . The uncertainty linked to our treatment of in- medium effects is reflected 
by the band corresponding to —250 < Ato£i(Tc) < —80 MeV, with stronger in-medium 
effects resulting in a smaller J /ip yield. 

5. Conclusions 

We have presented a model for charmonium production in heavy-ion collisions 
incorporating in-medium effects on open-charm states, as inferred from Lattice QCD, 
within a kinetic rate equation which allows to comprehensively treat suppression and 
regeneration mechanisms during the course of the collision. We have found that QGP 
formation manifests itself by J suppression at SPS energies and by J /ip regeneration 
at RHIC, where run-4 data are expected to give important insights. In-medium 
effects have so far proved to be essential to understand the centrality dependence 
of the ip' /ip ratio at SPS. Our predictions can also be tested by upcoming NA60 
data for In-In collisions at SPS. Complementary studies of charmonium transverse 
momentum distributions, as well as charmonium and bottomonium production at 
LHC, will provide further scrutiny of the proposed approach. 
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